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Toward understanding intermolecular interactions governing self-association of 
proteins, the present study investigated a model protein, myoglobin, using a 
small-angle X-ray scattering technique. It has been known that removal of the 
heme makes myoglobin aggregation-prone. The interparticle interferences of 
the holomyoglobin and the apomyoglobin were compared in terms of the 
structure factor. Analysis of the structure factor using a model potential of 
Derjaguin-Laudau-Verwey-Overbeek (DLVO) suggests that the intermole- 
cular interaction potential of apomyoglobin is more attractive than that of 
holomyoglobin at short range from the protein molecule. 
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1. Introduction 

Precipitation or aggregation of proteins is often problematic 
when handling protein solution for research and commercial 
use. Developing methods to inhibit undesired aggregation of 
proteins is one of the critical issues in the biopharmaceutical 
field (Shire et al, 2004; Yang et al, 2010). To our knowledge, in 
protein design studies, we often come across water-insoluble 
artificial proteins (Isogai et al, 2005; Imamura et al, 2012). 
In contrast, a mechanism which prevents an uncontrolled 
aggregation is likely to exist in naturally occurring proteins 
(Isogai, 2006; Monsellier & Chiti, 2007). Understanding 
physicochemical mechanisms governing self-association of 
proteins will be useful for future protein engineering. In the 
present study, we used a small-angle X-ray scattering (SAXS) 
technique to analyze protein's interparticle interference, 
which yields an intermolecular interaction of a protein in 
solution through use of the theory for liquids (Tardieu et al, 
1999). We employed myoglobin as a model protein. Removal 
of heme is regarded as a kind of simple protein engineering. 
Interestingly, it was reported that the apomyoglobin, in which 
the heme is removed, is prone to aggregate, while holo- 
myoglobin, in which the heme is binding, is highly water- 
soluble (Fandrich et al , 2001 , 2003). Crystallization of proteins 
would be controlled by protein's self-interaction. To date, a 
three-dimensional structure of apomyoglobin by X-ray crystal 
analysis has not been uploaded to the Protein Data Bank, 
while holomyoglobin was the first success in the X-ray crystal 



analysis of the protein structures (Kendrew et al, 1958). 
Unravelling how removal of heme changes the protein- 
protein interaction of myoglobin will be a step toward 
understanding key factors for the aggregation and the solu- 
bility. 

2. Materials and methods 

2.1. Sample preparation 

Holomyoglobin from equine skeletal muscle was purchased 
from Sigma. Apomyoglobin was prepared according to the 
following methods (Teale, 1959; Hapner et al, 1968). 
10 mg mP 1 aqueous solution of myoglobin was adjusted to pH 
~1.5 with HC1 at 297 K. An equal volume of ice-cooled 2- 
butanone was mixed with the solution, and then the organic 
phase containing the extracted heme was discarded. This 
procedure was repeated three times. The aqueous phase was 
dialyzed against sodium bicarbonate (0.05 g and then 
ultrapure water (Milli-Q, Millipore). The apomyoglobin 
solution was finally lyophilized. The solutions of holo- 
myoglobin and apomyoglobin for the measurements were 
prepared by dissolving the powders into the ultrapure water. 
The pH of the holomyoglobin and apomyoglobin solutions 
were 7.5 and 6.7, respectively. The concentrations of both the 
proteins were determined by using an £ 410nm of 1.6 x 10 5 for 
holomyoglobin and £280 nm of 1.43 x 10 4 for apomyoglobin 
(Goto et al, 1990). 
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2.2. Small-angle X-ray measurement 

SAXS was measured using a SAXSess camera (Anton Paar; 
Graz, Austria) with a sealed-tube X-ray source (PANalytical, 
The Netherlands), operating at 40 kV and 50 mA, and line 
collimation. X-ray intensities were recorded using an imaging 
plate (100 um x 100 um pixel size; Fuji Film, Japan), which 
was read using a Cyclone scanner (PerkinElmer, USA). The 
X-ray wavelength, X, was 1.542 A (Cu Ka); the camera length 
was 264.5 mm. The image data were integrated into one- 
dimensional scattering intensities using the program Image.! 
(developed by Wayne Rasband, National Institute of Health, 
Bethesda, MD, USA) with the macros in Utah SAXS Tools 
written by Professor David P. Goldenberg (University of Utah, 
USA). The scattering parameter q is defined as q = AitsmOIX, 
where 28 is the scattering angle of X-rays. The width of the 
integration area (detector slit width) was 10 mm. Protein 
solution and water were measured at 293 K in 1 mm-diameter 
glass capillaries for 1 h. 

2.3. Calculation of a structure factor from an intermolecular 
potential 

The SAXS intensity of a protein, I(q), is described by 

I(q) = ckP(q)S(q), (1) 

where c is the protein concentration, k is a constant, P(q) is a 
form factor (protein's self-scattering) and S(q) is a structure 
factor (protein's interparticle interference). The SAXS 
profiles of the dilute protein solutions (0.5 wt%), at the 
concentration of which the interparticle interferences are 
negligibly small, yields kP{q) by an indirect Fourier transform 
program GNOM (Svergun, 1992) involving correction of the 
slit smeared effect. S(q) can be connected to the inter- 
molecular potential of a protein molecule, V(r), by solving the 
Ornstein-Zernike equation with an appropriate closure. In 
the present analysis, we used the closure obtained from a 
random phase approximation, under which S(q) is described 
by 



S(q) = S 0 (q)[l + (n/k B T)S 0 (q)V(q)Y 



(2) 



(Hansen & McDonald, 1976; Narayanan & Liu, 2003), where 
So(q) is the structure factor of the reference system, n is the 
protein's number density, k B is the Boltzmann constant and 
V(q) is the Fourier transform of V(r). S Q (q) has been evaluated 
using the empty core model (Croxton, 1975; Kelkar et al, 
1992). The present study employed the Derjaguin-Laudau- 
Verwey-Overbeek (DLVO) model potential (Verwey & 
Overbeek, 1948), in which V(r) is expressed as 



V(r) = V m (r) + V c (r) + V AY (r), 



(3) 



where r is the distance between the protein molecules. The 
term Vus(r) is the hard sphere potential given by 



oo, 
0, 



r < a, 
r > a, 



(4) 



where a is the diameter of the protein. The terms V c {r) and 
^ ayW are the screened repulsive Coulomb potential and the 
attractive potential, respectively [here, V AY (r) is assumed to 
be a Yukawa-type potential], 



V c (f) = 



Z 2 e 2 exp[— k{t — a)] 



e(l + 0.5/ar) 2 



V A ^{r) = J{G/r) exp 



-{r-a) 



d 



(5) 



(6) 



where Z is the net charge on the protein, e is the elementary 
charge, e is the dielectric constant of the medium, k is the 
reciprocal Debye-Huckel screening length, / is the depth of 
attractive potential at r = a, and d is the range of the attractive 
potential. S(q) can be simulated, according to the equations 
described above and the given values for the parameters (Z, a, 
J and d). In the present study, we used Z = 1 and a = 41.2 A, 
the values of which were from the literature (Longeville et al., 
2003). The procedure for the present analysis was performed 
by using IGOR Pro version 6.22A (Wavemetrics, USA). 

3. Results and discussion 

Fig. 1 shows the observed SAXS profiles of holomyoglobin 
and apomyoglobin at 0.5 and 6.3 wt%, where the intensities 
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The normalized SAXS intensities of (a) holomyoglobin and (b) apomyoglobin at 6.3 wt% (crosses) and 0.5 wt% (open circles). The solid lines are 
theoretical fits oiI s (q) and the dashed lines are the calculated S(q), where the / and d values (0.7 fc B Tand 12 A for holomyoglobin, 14 fc B 7"and 0.7 A for 
apomyoglobin) were determined by the non-linear least-squares Levenberg-Marquardt algorithm. 



920 Hiroshi Imamura et a/. • Modulation of the intermolecular interaction of myoglobin 



/. Synchrotron Rad. (2013). 20, 919-922 



diffraction structural biology 




Maps of the residual sum of squares between the experimental I s (q) and the calculated I s (q) within ranges of / (0.01-9fc B 7) and d (1-15 A) for (a) 
holomyoglobin and (b) apomyoglobin. The darker portions in the map represent the values which efficiently reduce the residuals, (c) The protein- 
protein interaction potentials, V(r), are calculated from the same values of / and d in Fig. 1 (full line: holomyoglobin; dashed line: apomyoglobin). 



are normalized by the protein concentrations. For both the 
6.3 wt% solutions, the intensities at low scattering vectors (q < 
~0.1 A -1 ) were depressed compared with that for 0.5 wt% 
solutions, which are due to the interparticle interference, 
i.e. the structure factor. We analysed the structure factor 
component in the SAXS intensities based on the DLVO model 
interaction potential and the equations above. A slit-smeared 
SAXS intensity, I s (q), was simulated from I(q), which is from 
the determined kP(q) and the theoretically derived S(q) with 
equation (1), and the instrumental beam profile used as a 
weighting function (Glatter & Kratky, 1982). To find the 
values of the parameters (/ and d), which reproduce the 
experimental I s (q), we compared the experimental I s (q) with 
the calculated I s (q), varying the values of / (0.01-9 k B T) and d 
(1-15 A). The maps of the residual sum of squares are shown 
in Figs. 2(a) and 2(b). It is found that smaller / and larger d 
tend to reduce the residuals for holomyoglobin, while larger / 
and smaller d tend to reduce the residuals for apomyoglobin. 
The protein-protein interaction potentials, V(r), calculated 
with the optimal values of / and d, which give the best fits, are 
depicted in Fig. 2(c). At small r (r < ~45 A), the potential of 
apomyoglobin is more attractive than that of holomyoglobin, 
due to the larger / value. This could explain the irreversible 
aggregation-prone property of apomyoglobin (Fandrich et ah, 
2003). In addition, the attractive potential of holomyoglobin 
is relatively small negative at small r (r < ~45 A), but still 
remains at larger r (r < ~55 A) due to the larger d value. The 
attractive potential allows the proteins to be close to each 
other, but to be separated because the magnitude is compar- 
able with thermal energy (~k B T) at short range. This seems to 
underlie the highly water-soluble and less irreversible aggre- 
gation-prone properties of holomyoglobin. 

4. Summary 

In the present study we investigated the role for the heme in 
the self-interaction of myoglobin molecules by SAXS. From 
the fitting of the experimental SAXS data with the simulated 
SAXS intensities based on the DLVO model and the equa- 
tions in the liquid theory, it is suggested that change in the 



SAXS profile by removal of the heme could be explained 
by an increase in the intermolecular attractive interaction 
potential at short range from the protein molecule. We remark 
that the potential we can obtain depends on the selected 
model. Therefore, we progress a study developing a method 
to determine the intermolecular potential of proteins from 
experimental structure factors without assuming any model 
potential by applying an integral equation theory. 
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